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We study modulational instability of two-component Bose-Einstein condensates in an optical 
lattice, which is modelled as a coupled discrete nonlinear Schrodinger equation. The excitation 
spectrum and the modulational instability condition of the total system are presented analytically. 
In the long-wavelength limit, our results agree with the homogeneous two-component Bose-Einstein 
condensates case. The discreteness effects result in the appearance of the modulational instability 
for the condensates in miscible region. The numerical calculations confirm our analytical results 
and show that the interspecies coupling can transfer the instability from one component to another. 
PACS numbers: 05.45.-a, 03.75.Lm 



I. INTRODUCTION 

I The modulational instability (MI) is a general phe- 
nomenon in nonlinear wave equations and can occur in 

I various physical systems, such as fluid dynamics, plasma 
physics, and nonlinear optics. Due to the interplay be- 
tween nonlinearity and the dispersive effects, a weak per- 
turbation on a plane wave may induce an exponential 
growth, and result in the carrier wave break up into a 
train of localized waves [1]. Recently, the MI of Bose- 
Einstein condensates (BEG) has attracted much inter- 

I est, especially for the single-component BEG in an op- 
tical lattice [2,3]. In the superfluid regime, the BEG in 
an optical lattice can be modeled by a discrete nonlinear 
Schrodinger equation [4,5]. Following Ref. [6], the MI of 
a BEG in a deep optical lattice has been studied both 
theoretically [7-10] and experimentally [11,12]. 

Gompared with the single-component BEG system, bi- 
nary mixtures of BEGs have much richer structures [13]. 
In fact, the excitation spectrum of a homogenous two- 
component BEGs has been studied in Refs. [14-16]. It 
was shown that the MI depends strongly on the sign 
of ~ A1A2, where A12 is the interspecies interac- 
tion strength, Ai and A2 are the intraspecies scatter- 

' ing strengthes, respectively. In Ref. [17], the MI of the 
two-component BEGs in a deep optical lattice was stud- 
ied with a Josephson-like coupling term. However, their 
works were limited on a fixed perturbation wave number 
(i.e., q = tt). In addition, they did not compare their re- 
sults with that of the homogenous two-component BEGs 
case [14-16]. 

In this paper, we study the MI of two-component BEGs 
in an optical lattice. The explicit expression of the ex- 
citation spectrum is presented analytically. In order to 
outline the effects of the discreteness, we compare the MI 
condition of the modulated plane waves with that of the 
homogenous two-component BEGs case, in which the MI 
can occur only for A^2 > ^i-'^2, i-e., the BEGs in phase 
separation region [14-16]. Our results show that, the MI 
of the two component BEGs in an optical lattice depends 
not only on the sign of — A1A2, but also on the wave 



number k of the carrier waves. For A?2 < A1A2, the MI 
of the miscible condensates can also take place when the 
condensates have relatively large wave number k > 7r/2. 
The effect of the interspecies coupling is investigated by 
numerical calculations. Our results show that the cou- 
pling between the two condensates can transfer the in- 
stability from one component to another. 

II. THEORETICAL MODEL 

We consider a cloud of BEG atoms which have two 
internal states labeled by |1) and |2) in an optical lat- 
tice. When the heights of the interwell barriers are 
much higher than the chemical potentials, the conden- 
sate wave functions can be expressed as a sum of wave 
functions localized in each well of the periodic potential 
^'^(r, t) = Vn J2j V'j,o-(i)0<T(r - Tj), where N is the to- 
tal number of condensate atoms, and (/)cr(r — r^) is the 
spatial wave function localized in the jth site. With the 
help of the tight-binding approximation [4] , the probabil- 
ity amplitudes ipj a obey a set of coupled DNLS equations 

dip ■ (J 

+ {Aa,a\lpj,a\^ + Ka'\^j,a'\^) ^j,a, (1) 

where K^, and A^^i with cr = 1,2 describe the nearest- 
neighbor hopping term and the on-site atomic collisions, 
respectively. The stationary solutions of the DNLS equa- 
tions (1) are plane waves -iAj.o- = ■0^'^-'e**^'^^ with 
the chemical potentials of the two-component condensed 
atoms 

fi„ = cos(fc) + A,,„|vS°f + (2) 

Here, the chemical potentials fia sue to be determined 
by the normalization condition o- IV'i.o-P = 1- The 
terms —2KcrCos{k) in Eq. (2) are the kinetic energies 
aroused from tunnelling. It should be noted that, we 
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study the BECs system in the supcrfluid regime with 
{N/M)Kcr > A„^cr' ( M is the lattice number), and thus 
the dynamics can be well described by the DNLS equa- 
tions. 

The excitation spectrum can be obtained by the 
stability analysis of the plane waves. We introduce 
small perturbations to the carrier wave as: ipj.a- = 



phonon modes are taken in the traditional form Sijjj,^ = 
y^gi(gi-a;t) ^ y*g-i{qj-<^t) ^ ^-^j^ ^^le Constants and Va- 
to be determined. Then we obtain the excitation spectra 
as 



;W = sin(fc) sm{q) ± ^e<,(eg + A^), 



(3) 



where the small amplitude where 



A, = Ai(^i' 



(0)n2 



(0)n 



(4) 



In the above derivations, we have assumed ij^a to be 
real, and taken Ki = K2 = K, then eg,i = eq,2 = 
€q = 4K cos{k)sm^{q/2) [18]. In the absence of the in- 
terspecies interactions (A12 = 0), A„ = 2Ai{ip^''y or 

2A2(V'2°^)^- The two-component BECs are fully decou- 
pled, and the spectra are the same with single-component 
case [9,18]. According to Rcfs. [6,8], the plane waves of 
the condensates are modulationally unstable when the 
eigenfrequencies becomes imaginary, i.e., eq{eq+A„) < 0. 
In this case, the system undergoes exponential growth 
with the growth rate determined by the imaginary part 
of Uq^J. Prom Eq. (4), we know that: (i) if Afj < A1A2, 
i.e., the case of the two components being miscible, both 
Ai and A2 are positive; (ii) if Afj > A1A2, then Ai < 
and A2 > 0. Within the second scenario, the conden- 
sates tend to separate spatially, namely the phenomena 
of phase separation [19,20]. The sign of (Afj — A1A2) 
is of crucial important in discussing the stability of the 
eigenmodes. 



III. MODULATIONAL INSTABILITY 



cos(A;) < and + A^- > 0. For simplicity, we take 
= 1^2'^ = V'o, then Act = fi^^'oi where 

^)ct = Ai + A2 + (-l)"-\/(Ai - A2)2 + AAj^. (6) 

Therefore, the miscible BECs may exhibits the MI when 
7r/2 < k < 37r/2 and Vo > —^q/^a- More specially, 
both of the two-component with relative large wave num- 
ber {k > 7r/2) will be unstable for arbitrary q provided 
that tpQ > ipQ = 4K/ili. Below the critical am- 
plitude, the MI region depends on the perturbation's 
wave number q. In Fig. 1 (a), we consider the case 
tpQ < i/iQcri where the atomic collisions are taken as 
ai : ai2 :'a2 :: 1.007 : 1 : 1.01 [22]. We take Ai = A, 
then A12 = 0.993A and A2 = 1.00298A, which satisfy the 
miscible criterion. Our results show that there exist three 
different regions in the {q, k) plane. The dotted line in 
Fig. 1(a) separates the unstable regime into a fully un- 
stable and a partially unstable, in which one component 
is stable and the second one is unstable. 



In this section, we study the modulational instability 
by using Eq. (3) . To understand the effect of the discrete- 
ness clearly, we first study Eq. (3) in the long-wavelength 
limit [6]. Eq. (3) in the long- wavelength limit, i.e., fc — > 
and g — > 0, is reduced to: 

a;(±; = 2Kqk ± ^/K^{K^TA^, (5) 

which agrees with the usual expression obtained from 
continuous NLS equations [14-16] by replacing Kq^ with 
q^ /2m. It was shown that only immiscible BECs (with 
Af2 > A1A2) exhibit the MI [21]. 



A. The MI condition for the miscible BECs 

The MI condition for the two-component BECs are 
determined by the condition: eg(eg -|- A^) < 0. For 
the miscible BECs with Af2 < A1A2, the MI occurs for 



B. The MI condition for the BECs in phcise 
separation region 



To study the MI of the immiscible BECs, wc take the 
atomic collisions arc taken as ai : ai2 : 0,2 1-03 : 1 : 0.97 
[23]. Thus wc have A12 = 0.9709A and A2 = 0.9417A, 
which meet the condition Afj > A1A2. Unlike the previ- 
ous miscible case, the MI condition considered here de- 
pends on the wave number k. We find that for cos A; > 
(i.e., k < 7r/2 or A; > 37r/2), the MI condition is 
V'o > However, for 7r/2 < k < 3tt/2, it becomes 

tpQ > —eq/^2- Since <C ^2, the immiscible BECs 
will be unstable in the whole (q, k) plane provided that 
tpQ > ipQ = 4:K/\ili\. Below the critical amplitude, the 
MI region in the (g, k) plane is shown in Fig. 1 (b). The 
MI can take place for any wave number k. 
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FIG. 1. Region of the MI in the (q, k) plane for (a) miscible 
BEG Ai2 = 0.993A and A2 = 1.00298A, (b) immiscible BEG 
A12 = 0.9709A and A2 = 0.9417A. Other parameters are 
taken as JiT = 1, Ai = A = 100, and Vo = 1/(2M + 1) < Vo.cr- 

IV. NUMERICAL ANALYSIS 

The linear-stability analysis can determine the onset of 
MI, however it does not yield any dynamical information 
beyond the instability point. Therefore, we perform nu- 
merical simulations of the DNLS equations. The initial 
conditions are two modulated plane waves 

= V',-,2(0) = [A + a cos(gj)]e''^ , (7) 

in accordance with our previous discussions. The modu- 
lation amplitudes a are assumed small compared with the 
backgroimd amplitudes a = 0.05^. We consider the op- 
tical lattice with the total number of the sites M = 400. 
Prom the periodic boundary condition, the wave num- 
bers take the form of fc = 27r//M and q = 2tts/M, where 
we choose the unit lattice constant. 

We first consider the miscible BEC case. We take the 
amplitude as A = 1/^2M + 1 to insure the normaliza- 
tion condition. The wave numbers are chosen as {I = 50, 
s = 100), {I = 150, s = 100), and {I = 150, s = 50), 
which correspond to the points labeled by 1, 2, 3 in Fig. 
1(a). The first two points are stable, and the last one is 
in an unstable regime with the growth rate being 0.1863. 
Density of the first component \tpj,i{t)\^ is plotted in Fig. 
2. As expected from the analytical prediction, for the 
first point (fc = 7r/4, q = 7r/2), the modulated wave is sta- 
ble. However, for the second point (fc = 37r/4, q = 7r/2) 
as shown in Fig. 2(b), the density increase sharply at 



the final stage. It is not predicted by the linear stability 
analysis. The reason for this phenomenon is the com- 
plex interactions between three fundamental eigenmodes 
of wave numbers fc, k — q, and k + q contained in the 
initial state, and additional components of other wave 
numbers [6]. Fig. 2(c) shows the case that iVq^i is real 
and u}q^2 is imaginary, which means that the component 
1 is stable and the component 2 is unstable. However, 
the instability is transferred from one component to an- 
other due to the effects interspecies coupling, which in 
turn leads to the MI of total system. The appearance of 
spatially localized modes with large amplitude indicates 
the modulational instability of the miscible BECs. 




site index j 

FIG. 2. Density of component 1 for the miscible case 
(A12 < A1A2). Gomponent 2 is similar with component 1. 
(a) k = n/4, q = 7r/2, (b) fc = Sn/A, q = n/2, and (c) 
fc = 877/4, g = 7r/4. Other parameters are the same as Fig. 
1(a). 

In Ref. [17], Rapti et al. calculated the spatiotempo- 
ral dynamics of the two-component BECs with purely 
nonlinear coupling, where the wave number fc is taken as 
'7r/2. According to the linear-stability analysis, however, 
within this case both Uq^i and u)q^2 are real and equal 
to 2Ksm{q), which means that there is no the MI. The 
appearance of the large-amplitude excitations in their nu- 
merical simulation may originate from the complex inter- 
actions between the eigenmodes involved. 

We further study the MI of the two-component BECs 
for the immiscible case. As an example we calculate two 
points: {I = 150, s = 5) and {I = 150, s = 10) in the 
unstable regime. Both points have real cuq^i and imagi- 
nary Wq,2- The MI growth rates are 0.0458 and 0.0902, 
respectively. The instability growth time for the second 
point is twice to that of the first point with the second 
one. Our numerical results shown in Fig. 3 confirm the 
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analytical prediction. Other simulations in the unsta- 
ble regime show that the increase of q in the monotonic 
regime (0 < g < tt) can enhance the growth rate, i.e., the 
instability growth time becomes more and more shorter 
with the increase of q. 




site index j 



FIG. 3. Density of component 1 for the immiscible case 
(Af2 > A1A2). Component 2 is similar with component 1. 
(a) A; = 37r/4, g = 7r/40, (b) k = 37r/4, q = 7r/20. Other 
parameters are the same as Fig. 1(b). 

V. CONCLUSIONS 

We have studied the modulational instability of two- 
component BEG in the superfluid regime. The explicit 
expression of the excitation spectrum is presented ana- 
lytic;ally. In the long- wavelength limit, our results can 
recover previous results of two-species BECs in homoge- 
neous case. From the Bogoliubov dispersion relation we 
studied the MI condition of the modulated plane waves. 
The discreteness effect leads to the appearance of the MI 
for the two-component BECs not only in the miscible 
region but also in the phase separation region. 

For the miscible condensates, the MI can occur for rel- 
atively large wave number of the condensates k > 7r/2. 
For the smaller tpQ compared with the critical density 
TAo.crj the MI condition depends also on the perturbation 
wave number q. The MI conditions in the two-component 
BECs with the phase separation can take place for any 
k. However, unlike to the miscible case, only one of the 
species can be modulationally unstable for any given q 
and k. To confirm the analytical results, we also per- 
formed numerical calculations. We find that the inter- 
species coupling can transfer the instabilities from one 
component to the other, which lead to the MI of the 
whole system. 



In summary, with the help of linear-stability analy- 
sis, we have studied the stability of the modulated plane 
waves, which is essential to predict the existence of non- 
linear localized modes, such as dark solitons, in the dis- 
crete nonlinear Schrodinger equation [24,25]. 
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